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ABSTRACT 



We report the discovery and the Rossiter-McLaughlin effect of Kepler-8b, a transit- 
ing planet identified by the NASA Kepler Mission. Kepler photometry and Keck-HIRES 
radial velocities yield the radius and mass of the planet around this F8IV subgiant host 
star. The planet has a radius Rp = 1.419 Rj and a mass, Mp = 0.60 Mj, yielding 
a density of 0.26 gem -3 , among the lowest density planets known. The orbital pe- 
riod is P = 3.523 days and orbital semimajor axis is 0.0483lg;gg^AU. The star has 
a large rotational vsini of 10.5 ± 0.7 kms _1 and is relatively faint (V ~ 13.89 mag), 
both properties deleterious to precise Doppler measurements. The velocities are indeed 
noisy, with scatter of 30 ms _1 , but exhibit a period and phase consistent with the 
planet implied by the photometry. We securely detect the Rossiter-McLaughlin effect, 
confirming the planet's existence and establishing its orbit as prograde. We measure 
an inclination between the projected planetary orbital axis and the projected stellar 
rotation axis of A = —26.9 ± 4.6°, indicating a moderate inclination of the planetary 
orbit. Rossiter-McLaughlin measurements of a large sample of transiting planets from 
Kepler will provide a statistically robust measure of the true distribution of spin-orbit 
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orientations for hot jupiters in general. 

Subject headings: planetary systems — stars: fundamental parameters — stars: indi- 
vidual (Kepler-8, KIC 6922244, 2MASS 18450914+4227038) 



1. Introduction 



To date, 90 "hot j 
around Sun-like stars 



upiters" — gas gian t planets with periods <10 days — have been detected 
Torres et al.l 120081 ) . The front-running formation scenario supposes that 



these planets did not form where they reside today, close the host star, because the inner region s 
of protoplanetary disks have inadequate surface densities and high temperatures (|Lin et al.l ll996). 
Instead hot jupiters are presumed to form several astronomical units (AU) from their host stars 
followed by subsequent migration inward to their current locations. One likely migration sce- 



nario involves tidal intera ctions between the planet and a remaining gaseous disk (ILin et al. 



1996 



Moorhead &: Adamsll2008l ) , causing the planet to spiral inward while maintaining its nearly circular 
orbit that is co-planar with the disk. Alternati vely, migration may o ccur by N-body gravita tional 



interactions, such as planet-planet scattering (IRasio fc Ford 1996; 



namical relaxation (jPapaloizou fe Terquem 



2001 



Adams fc Lau ghlin 



with a distant object , and damped later by tidal friction (jHolman et al .11 199 71 : iFabrvckv h Tremaine 



2007 



Chatteriee et al. Il2008h . dy- 
2003|) . and Kozai interactions 



Wu et alj|2007l : iNagasawa et al.ll2008l ). Thus, measurements of both the orbital eccentricities 
and the orbital inclinations relative to the star's equator offer diagnostics of the original migration 
process. 

We assume that planets form in protoplanetary disks with the stellar spin and planetary 
orbital axes aligned. If so , the nearly adiabati c tidal interactions between planets and disks would 
maintain the alignment (IWard &: Hahnl Il994l ). In contrast , few-body gravitati o nal in teractions 
would typically cause misalignments. Few-body models by lAdams Laughlin I (|2003l ) predict a 
final inclination distribution for dynamically relaxed planetary systems that peaks near 20° and 
which extends to inclinations as high as 85°. Kozai interactions between a planet and an outer 
body (star or planet) result in a wide distribution of final orbital inclinations for the inner pla net, 
including retrograde orbits (IFabrycky <fe Tremaind 120071 : IWu et al.1 120071 : INagasawa et al.ll2008l ). 



The Rossiter-McLaughlin (R-M) effect offers a way to assess quantitatively the spin-orbit align- 
ment of a planetary system by measuring the Doppler effect of the star's light during a planetary 
transit. As the planet blocks receding portion of a rotating star's surface, the spectrum from the un- 
obscured surface has a net Doppler shift toward shorter wavelengths, and vice versa fo r blocking the 
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About 2/3 of the 18 planetary systems measured by the R-M effect have an orbital plane well 
aligned with the star's equatorial plane, as projected onto the sky, giving A near 0°. This alignment 
is as expected from simple migration theory due to gentle loss of orbital energy to the gas in the 
protoplnaetary disk (|Lin et al.l ll996). H owever, six exoplanetary systems show a significant spin- 



orbit misalig nment, namel y HP 806 06 (W inn et al 



Winn et al 



Gillon 2010), WASP-14b (Johnson et al 



2009), HAT-P-7b (Winn et al. 



Wasp- 17b (| Anderson et al. Il2009h . 



200 



* 



2009 



2009 



Joshi et al 



Narita et al. 



Moutou et al. 2009: Pont et al. 



2009), XO-3b (IHebrard et al 



2009al ). CoRoT-1 (jPont et al. 



2009: 



2009 



2009) 



The variety of alignments support the bimodal distribution found bv lFabrvckv Winnl (120091 ) . 
Interestingly, all four misaligned systems contain quite massive planets, above 1 Mj. This cor- 
relation may be re lated to the association of massive planets with higher orbital eccentricity 
(jWright et al,ll2009l ). as both eccentricity and inclination may arise from perturbations of plan- 
ets from their original circular orbits. But there are two mas sive planets on eccentric orbits for 



which A appears to be con sistent with zer o, namely HP 17156b (ICochran et al 



2009 



Narita et al1l2009bh . and HATp-2b dwinn et al.1120071 : lLoeillet et al.lhoosh . There is currently 



2008; 



Barbieri et al. 



no dominant and secure explanation for the misaligned or eccentric hot jupiters. 

Here, we present the first detection of the Rossiter-McLaughlin effect from a planet detected by 
the Kepler Mission. As this mission is expected to detect dozens of transiting hot jupiters, Kepler 
offers an opportunity to provide a statistically robust measure of the distribution of spin-orbit 
angles, and to correlate that angle with other physical properties of the systems. 



2. Kepler Photometry 



Nearly continuous photometry in a 100 square degrees field near Cygnus and Ly ra was carried 



out during 42 days by the Kepler spaceborne telescope, as described previously ( Borucki et al 



2010 



Borucki et al. 



2010 



Koch et al. 



2010: 



Jenkins et al 



2010; 



Batalha et al. 



2010; 



Gautier et al 



20ld ). The star Kepler-8 (= KIC 6922244, a = 18 h 45 m 09M5, 5 = +42°27'03 // .9, J2000, KIC 



r = 13.511 mag) exhibits a repeated dimming of 9.82 ± 0.22millimag, obvious against uncertainties 
in each 30 minute integration of 0.1 millimag. The light curve for Kepler-8 is plotted in Figure 
[TJ The numerical data are available ele ctronically in the online edition of the journal. A modest 
amount of detrending has been applied (jKoch et al.ll2010i : iRowe et al.ll2010l ) to the time series. 



We detect no systematic difference between alternating transit events at 50 /imag levels, ruling 
out nearly equal components of an eclipsing binary star, (see Fig.[TJ). We also see no evidence of 
dimming at the expected times of a secondary eclipse, which would be visible for most eclipsing 
binary systems of unequal surface brightness. The photocenter shows no displacement astrometri- 
cally above millipixel levels (0.5 mas) during times in and out of transit as would be seen if there 
were a background eclipsing binary masquerading as a transiting planet. Thus the photometric 
and astrometric non-detections of an eclipsing binary support the planet interpretation for the re- 
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peated transit signatures. Moreover, the shape of the photometric transit is adequately fit with a 
planet-transit model further supporting the planet interpretation. 

We fit the light curves by solving for a/Rstan the density of the star, and the ratio of planet 
to stellar ra dius. We follow the method for es t imati n g stellar rad i i and othe r stellar parameters 



descri bed by lSozzetti et al.1 (120071 ) . Bakos et all (120071 ). IWinn et al.1 (120071 ) and lCharbonneau et al 
(|2007n . This method extracts physical properties directly from the light curve, geometry, and 



Newtonian physics, and it uses the greater orthogonality of parameters to yield more robust fits to 
observables (jTorres et al.ll2008l ) . We fit explicitly for a/R^, b , and r/it^using a pro cedure developed 
by one of us (J.Rowe) and described bv lKoch et al.l (|2010l ); iBorucki et al. I (|2010i ). 



W e begin with an LTE spectroscopic analysis (IValenti Piskunovi 1 19961 : IValenti Fischer 
2005) of a high resolution template spectrum from Keck- HIRES of Kepler-8 to derive an effective 
temperature, T e ff= 6213 ± 150 K, surface gravity, log g= 4.28 ± 0.10(cgs), metallicity, and the 
associated error distribution for each of them. The multitude of Yale- Yonsei stellar evolution models 



(jDemarque et al .1120041 ) are constrained by both those LTE measurements and by the stellar density 
that stems directly fr om the orbital period, the fractional dim ming during transit, and measures 



of transit durations (jSozzetti et al.l 120071 : iBrown et al.l 120101 ) . By Monte Carlo analysis, those 



photometric and spectroscopic constraints and their uncertainties establish the probabiliby density 
contours among the evolutionary tracks where the star may reside. We iterate the self-consistent 
fitting of light curves, radial velocities, and evolutionary models until a domain of stellar mass, 
radius, and age is identified. That domain encompasses a range of evolutionary states that satisfy 
all of the constraints within their error distributions. The resulting mass and radius of Kepler-8 
are given in Table 1, along with other associated stellar properties such as luminosity and age. 



3. A Background Eclipsing Binary: Follow-up Imaging 

We carried out extensive tests of the possibility that the apparent photometric transit was 
actually caused by a background eclipsing binary star within the photometric aperture of radius 
~8 ". We obtained images with 0.8 arcsec seeing with the Keck telescope HIRES guider camera and 
the bg38 filter to search for stellar companions that might be eclipsing. This filter combined with 
the CCD detector have a response similar to the V plus R bands, in turn similar to the bandpass of 
the Kepler photometer. This Keck image is shown in[2l There is one star having 0.0075 the flux of 
the main star (V and R band) that resides 3.8"northwest of Kepler-8. This background star resides 
within the Kepler aperture and could conceiveably be an eclipsing binary that is masquerading as 
a transiting planet. 

But this faint neighboring star cannot be the source of the photometric transit signature for 
two reasons. If the background star were the cause of the observed 1% dimming, the photocenter 
centroid shifts woul d be 3-20 millipixel on the Kepler CCD. Instead, astrometric measurements in 



and out of transit (IJenkins et al.ll2010i ) show shifts of no more than 0.1 millipixel in and out of 
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transit. Figure [3] shows both the flux and astrometric photocenter (centroid) of the Kepler images 
during quarter 1 month of data. We applied a high-pass filter to remove non-transit signatures 
on timescales longer than 2 days. At times of transits there are no displacements in either the 
row or column direction at a level above a millipixel. This indicates that any background eclipsing 
binary would need to be well within ~0.1 pixels or 0.4"of the target star in order to explain the 
photometric transit signals. 

To hunt further for background eclipsing binaries we plotted flux versus the photocenter cen- 
troids in both the row and column directions, as shown in Figure HI These so-called "rain plots" 
would reveal an eclipsing binary as a "breeze" in the centroids to the left or right as the flux drops. 
No such breeze is detected at a level near 0.1 millipixel, ruling out all but eclipsing binaries located 
within a few tenths of an arcsec of Kepler-8. Finally, we looked directly at the Kepler images taken 
both during and out of transit to detect motion of the centroid of light, as shown in Figure 03 
We formed the difference of the images in and out of transit to detect astrometric displacements 
associated with the flux dimming, as would occur if a neighboring eclipsing binary were the cause. 
Those difference images show no shift of the photocenter. We conclude that the transit photometry 
with its 1% dimming cannot be explained by any eclisping binary companions beyond 1 arcsecond 
of Kepler-8. 

Furthermore, the flux ratio of the two stars is only 0.0075, making it impossible for the back- 
ground star to cause the 1% photometric dimming. Even if that background star were to vanish, 
the total flux would decline by less than the observed 1%. We conclude that the 1% photometric 
dimming is not caused by an eclipsing binary star within the Kepler photometric aperture, from 1 
- 10 arcsec of the target star. 

To hunt for additional stars in the field within an arcsec of Kepler-8 we used both speckle 
imaging and adaptive optics imaging. Kepler-8 was observed at the Palomar Hale 20 0-inch telescope 



on 09 Sep 2009 UT with Palomar near-infrared adaptive optics system [PHARO] (jHavward et al 



2001 



). The PHARO instrument was utilized in the J-band filter with the 25 mas per pixel (25"FOV) 
mode. The source was observed with a 5-point dither and an integration time of 2.8 seconds 
per frame. The dither pattern was repeated 7 times for a total on-source integration time of 98 
seconds. The average uncorrected seeing during the observations at J-band was 0.65", and the 
average AO corrected images produced point spread functions that were 0.09"FWHM. There are 
two faint sources within 4"of Kepler-8; these objects are 7 and 8.4 magnitudes fainter than the 
primary target and are too faint to produce the ifep/er-observed transit and centroid shift. The 
AO imaging detected no sources at J-band down to within 0.1" — 0.2" of the primary target that 
are within Am ~ 6 — 7 mag at J. 

Speckle observations of Kepler-8 were made on 1 October 2009 (UT) at the WIYN observatory 
located on Kitt Peak. The observations were made using the WIYN speckle camera during a night 
of very good seeing (0.62 ") and under clear conditions. We used a narrow bandpass 40 nm wide 
centered at 692 nm. The Kepler speckle program obtains observations of both double and single 
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standard stars throughout the night. We use a robust background estimator on the reconstructed 
images to set a limit for the level of companion star we should detect if present. The speckle 
observations show that Kepler-8 has no companion star between 0.05 and 2.0 arcsec within a delta 
magnitude of <4.2 mags. These AO and speckle observations effectively rule out the possibility of 
an eclipsing binary star to within as close as 0.1 arcsec from Kepler-8. 



3.1. Covariance of Inclination, Limb Darkening, and Stellar Parameters 



As Kepler-8 has a large impact parameter, the solution is quite sensitive to errors in limb- 
darkening. The derived impact parameter and inclination angle are directly related to the assumed 
limb-darkening law, poorly known for the wide Kepler bandp ass. We adopted an adhoc approach 
to estimate the limb-darkening parameters (jRowe et al.ll2010l ) as follows. We noted that fits to the 
observed Rossiter-McLaughlin effect (see section 6), coupled with the measured vsini= 10.5 ± 0.7 
km s _1 offered a constraint on the inclination angle and the impact parameter, and hence on the 
inf erred stellar radiu s, for a given limb-darkening law. To solve the genereal case for eccentric orbits 



sec 



Pal et al. I (|2010f ). 



For the cases of Kepler-4 and K epler-5 dBorucki et al. Il2010l : iKoch et al.ll2010l ). we discovered 
that a model using tabular values (jPrsa Zwitter 1 120061 ) for limb darkening over- predicted the 
amount of curvature in the variation of flux as a function of time during transit. We therefore 
modified the limb-darkening parameters by fitting the three known exoplanets, TRES-2, HAT-P-7, 
and HAT-P-11, and using published values for their stellar and planetary parameters. We linearly 
interpolated over the values of effective temperature to derive a superior measure of the true limb- 
darkening law for the Kepler bandpass. The photometry from new transiting planets found by 
Kepler, along with IR photometry that is less sensitive to limb darkening, will allow us to refine 
the limb-darkening law for the Kepler bandpass in the coming months. The resulting fits to the 
photometri c transits of Kep ler-4, 5, 6, 7, and 8 were all improved with our new limb-darkening 
treatment ([Rowe et al.ll2010l ). 



For Kepler-8, the resulting impact parameter and inclination were 6=0.724 ± 0.020 and i = 
84.07±+0.33 both about 10% smaller than we had obtained by using the first-guess limb-darkening 
law. The transit duration in turn dictates the best-fit stellar radius, and after iterating with the 
Yale-Yonsei models, the resulting stellar radius is R = 1.486^q'qq2 ^0- The stellar mass has a value 
within a narrow range, M = 1.213^Q Qg3 Mq, again derived from iteration between the fit to the 
light curve and to the Yale-Yonsei models. Additional information on stellar mass comes from the 
projected rapid rotation of the star, vsini= 10.5 ± 0.7 kms" 1 . This v sin zis faster than typical 
rotation for stars at the lower end of the mass range quoted above, especially as some evolution and 
associated increase in moment of inertia has occurred. Thus we marginally favor the upper half of 
the quoted mass range, i.e. M > 1.213 M . The stellar gravity is logg= 4.174 ± 0.026(cgs), and 
the approximate age is 3.84 ± 1.5Gyr, both indicating a star nearing the end of its main sequence 
lifetime. 
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4. Radial Velocities 

We took hig h resolution spectra of Kepler-8 using HIRES on the Keck I 10-m telescope 



(jVogt et al.l Il994j ) . We set up the HIRES spectrometer in the same manner that has been used 



consistently for 10 years with the California planet search (jMarcv et al. 1120081 ). We employed the 
red cross-disperser and used the I2 absorption cell to measure the instrumental profile and the 
wavelength scale. The slit width was set at Of! 87 by the "B5 decker", giving a resolving power 
of about 60,000 at 5500ABetween 2009 June 1 and 2009 October 31 we gathered 15 spectra of 
Kepler-8 out of transit (See Figured]). Typical exposure times were between 10 - 45 min, yielding 
signal-to-noise of 20 - 40 pixel -1 . These are quite low S/N ratios with which to attempt precise 
Doppler measurements. We tested the Doppler precision from such low SNR spectra on standard 
stars, notably HD 182488 and HD 9407, and found that the expected errors are ~5 ms , as ex- 
pected from Poisson statistics of the photons, with no systematic errors above 1 ms -1 , as usual for 
the iodine technique. On 2009 October 29 (UT) we obtained nine spectra of Kepler-8 during transit 
(the last taken during egress). The observations were taken between airmass 1.9 and 3.8, while the 
star was setting, the last exposure occurring at hour angle 5 hr 35 min with extreme atmospheric 
refraction and dispersion. Again, tests with standard stars observed at such high airmass show no 
systematic errors in the Doppler measurements above 1 ms -1 . 

We carried out careful reduction of the raw images, including cosmic ray elimination and 
optimal extraction of the spectra, to minimize the backg round moonlight. We performed the 



Doppler analysis with the algorithm of lJohnson et al.l (|2009l ) . We estimated the measurement error 



in the Doppler shift derived from a given spectrum based on the weighted standard deviation of the 
mean among the solutions for individual 2 A spectral segments. The typical internal measurement 
error was 7-10 m s . The resulting velocities are given in Table 2 and plotted in FigureEl 

The actual uncertainties in the velocities are certainly closer to 20-25 ms -1 for several reasons. 
The spectra have such low photon levels, ~200 photons per pixel in the raw CCD images, that 
cosmic rays and background sky are significant noise sources. The latter was especially problematic 
as we used an entrance slit with dimensions 0.87 x 3.5 arcsec, not long enough to separate the 
wings of the star's point spread function from the background sky. Nearly all HIRES spectra were 
taken with the moon gibbous or full, and about half of the nights had moderate cirrus that scatters 
moonlight into the slit. A few observations made with a 14 arcsec slit revealed that 1-3% of the light 
came from the moonlit sky in typical observations. Simulations with stellar spectra contaminated 
by moonlit sky suggested that errors of ~10 ms -1 would accrue, no doubt systematic as well. Errors 
in the velocities are large because of relatively rapid rotation, v sm.i= 10.5 ± 0.7, broadening the 
absorption lines by 5x that of the slowest rotating stars. Doppler errors increase linearly with 
line widths. Thus, while the S/N=40 might be expected to yield velocities of ~5 ms -1 , the high 
■usini increases that error to ~25 ms -1 . Indeed, the best orbital fit to the velocities, shown in 
Figure[6l exhibit discrepancies of 20-30 ms -1 . We estimate that the true errors in velocities are 
thus ~25ms _1 . We have accounted for these errors by adding a "jitter" of 25ms _1 to the internal 
velocity errors. Those augmented uncertainties are reflected in Table 2. 
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We carried out a Levenberg-Marquardt least-squares fit of a Keplerian, single-planet model to 
the observed velocities. In all models the orbital period and the time of mid-transit was constrained 
to be that found in the photometric fit. We first assumed a circular orbit leaving only two free 
parameters in our fit, namely the velocity amplitude and the gamma velocity of the system. The 
best-fit model is overplotted in Figured The best-fit value of the amplitude is K = 68.4±12.0 m s _1 . 
This, coupled with the adopted stellar mass of 1.2131q oes3 ^0> yields a planet mass 

We also carried out fits in which the eccentricity was allowed to float. The best-fit model 
has an eccentricity = 0.24, K = 71.4 ms" 1 , yielding a planet mass of 0.62 Mj. The RMS of the 
residuals is 39 ms _1 compared with 40 ms -1 for the circular orbit, rendering the eccentric orbit no 
better than the circular. Indeed a bootstrap Monte Carlo estimate of the parameter uncertainties 
gives a formal error in eccentricity of 0.16 . We conclude that the eccentricity is consistent with 
zero, but could be as high as 0.4. This large uncertainty stems from the high v sini of the star and 
it's relative faintness, V = 13.9 mag. 

A line bisector analysis showed no variation at a level of a few meters per second and no 
correlation with measured radial velocities. Thus, the velocity variation appears to represent actual 
acceleration of the center of mass of the star. 



5. Properties of the Planet Kepler-8b 

With the planet interpretation highly likely, the properties of the host star and the depth of 
transit directly yield a planet radius of 1.419^q Q5gi?j and a mass of 0.603t{j;JjjMj • The errors 
represent the 68%-probability domain of integrated uncertainty from all input measurements and 
models. The planet density is 0.261 ± 0.071 g cm _3 apparently placing it among the low density 
exoplanets. Assuming a bond albedo of 0.1, its orbital radius of 0.0483lQQQ ( j ) 2 AU, and the best-fit 
stellar luminosity from the Yale-Yonsei models gives an equilibrium temperature for the planet of 
1764 db 200 K, assuming rapid and complete redistribution of thermal energy around the planet's 
surface. This planet apparently is a member of the population of low density, bloated hot jupiters 
having thermal histories yet to be firmly understood. The planet properties are listed in Table 1. 



6. Rossiter-McLaughlin Effect 



The radial velocities ob tained during the tran sit of 2009 October 29 (UT) were modeled using 
the techniques described by lCochran et al.1 (|2008l ). The models adopted the system parameters in 
TableHJ an d searched for the value of A that best matched the observed velocities. In particular, we 
assumed the zero eccentricity orbit, as the K epler photometric lightcurve gives no indication of non- 
zero eccentricity. We used a IClaretl (|2004l ) four parameter limb-darkening law, although running 
similar models with both a linear and a quadratic limb-darkening law gave essentially the same 
results. We derive a value of A = —26.9 ± 4.6° for the angle between the projected stellar rotation 
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axis and the projected planetary orbital axis. The model fit to the data gives a reduced chi-squared 
Xq = 0.66 for the velocities obtained during the transit, indicating that the uncertainties in these 
velocities are probably slightly overestimated. 

The radial velocities measured during the transit and the Rossiter-McLaughlin effect model 
fit are shown in FigureO The observed asymmetric Rossiter-McLaughlin effect, with a positive 
deviation of larger amplitude and longer duration than the negative amplitude is a result of the 
combination of the large impact parameter and the non-zero value of A. A central transit would 
give a symmetric Rossiter-McLaughlin deviation, with the amplitude of the effect then depending 
on A. However, for an off-center tr ansit, the overal l shap e depends critically on both the impact 



parameter and on A (cf. Figure 2 of lGaudi fc Winnl (|2007l )). 



This Rossiter-McLaughlin analysis provides independent verification of the major results of 
the transit lightcurve analysis. The lightcurve analysis alone suggested a transit chord that is 
significantly off-center, with impact parameter, b = 0.724±0.020 and inclination, i = 84.07=L+0.33. 
This large impact parameter coupled with the non-zero value of A results in the planet transit 
blocking mostly the approaching (blue-shifted) half of the rotating star, as indicated in [8J The 
planet crosses to the receding (red-shifted) portion of the stellar disk just before the end of the 
transit. This causes the observed asymmetry in the Rossiter McLaughlin velocity perturbation 
during the transit. The observed amplitude of the R-M effect is in excellent agreement with 
the photometrically determined stellar radius and impact paramter, and with the spectroscopically 
determined v sini. The duration of the photometric transit agrees with the duration of the observed 
R-M velocity perturbation. This consistency between the properties of the system as derived from 
the transit photometry and from the R-M velocity perturbation gives a confirmation that both 
phenomena were caused by an orbiting, planet-sized companion to the star. There is no other 
explanation of the observed R-M variations that is also consistent with both the photometric and 
the radial velocity observations. 

This observation of the Rossiter-McLaughlin effect in Kepler-8 is important because it confirms 
that the observed lightcurve and radial velocity variations in Kepler-8 were indeed caused by an 
orbiting planetary companion. The amount of scatter in the radial velocity data giving the solution 
shown in Figure [6] is uncomfortably large, even given the lack of line bisector variations in these 
spectra. However, the R-M data absolutely confirm the existence of the planet. There is no other 
explanation of the observed R-M variations that is also consistent with both the photometric, 
astrometric, imaging, and the radial velocity observations. 



7. DISCUSSION 

We have carried out a wide array of observations of Kepler-8 that pinpoint the existence and 
properties of the exoplanet. The Kepler photometric measurements, with ~0.1 mmag errors in 30 
min intervals, are tightly fit by a model of a transiting exoplanet (Figure [1]). This good fit to the 
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inflections of the photometric data provide immediate support for the planet interpretation, with 
few plausible alternative interpretations except for a background, diluted eclipsing binary system 
having just the right brightness and radius ratio to mimic a planet. We note that, unlike ground- 
based transit work that has lower photometric precision, the Kepler photometry is so precise that 
blends of eclipsing binaries are more readily identified from the photometry alone. 

Nonetheless, to test the unlikely possibility of a blend, we carried out a battery of astrometric 
tests to hunt for an eclipsing binary. The resulting steadiness of the position of the primary star, 
Kepler-8, during and out of transit, argues against any eclipsing binary in the photometric aperture. 
We further carried out both adaptive optics imaging and speckle interferometric measurments to 
hunt for faint eclipsing binaries located within an arcsec of the Kepler-8. None was found, further 
diminishing the chance of such a masquerade. We followed with high resolution spectroscopy at 
both low and high signal-to-noise ratio, finding no evidence of double lines nor rapid rotation. 
Further support for the planet interpretation came from the precise radial velocities that varied in 
phase with, and had the same period as, the photometric light curve, further supporting the planet 
model and constraining the planet mass. 

Finally, the Rossiter-McLaughlin effect confirmed independently the planet interpretation and 
provided further geometrical information about the orbit, notably A = —26.9 ± 4.6°. A sketch 
of the star and planet's orbit are shown in (Figure [H]). Remaining unknown is the inclination of 
the star's rotation axis, as indicated in [HJ But continued photometry during the Kepler mission 
lifetime may reveal a photometric periodicity caused by the rotation of spots around the star. The 
resulting rotation period, coupled with the measured rotational vsini= 10.5 ± 0.7of the star, will 
allow the star's inclination to be measured, putting the R-M geometry on firmer ground. 

The suite of ground-based observations obtained for Kepler-8described above over-constrained 
a set of parameters that were obligated to mutually agree, including such subtle issues as the limb- 
darkening, the stellar rotation rate, and stellar radius, mass, and density. Thus, the R-M effect can, 
in general, be used to confirm the existence of a planetary companion in the case where standard 
follow-up observing procedures may provide somewhat ambiguous or equivocal results. Examples 
would be Kepler stars that are sufficiently saturated that measurement of photocenter shifts during 
transits is problematic, or stars like Kepler-8with large v sinithat make orbit determination via 
precise RV measurement difficult. 

Observations of the Rossiter-McLaughlin effect in dozens of transiting planets, discovered from 
the ground and from Kepler, offer an excellent opportunity to determine the distribution of orbital 
geometries of the short-period planetary systems in general. These planetary systems will continue 
to be studied in an extraordinarily unform and consistent manner with the same set of tools. The 
sample selection effects are extremely well understood and documented. The Kepler lightcurves 
are of unprecidented precision, allow sensitive searches for additional planets in the system via 
transit timing variations, and significantly tighter limits to be placed in the future on the orbital 
eccentricity. Coupling these data with R-M measurements of the orbital alignment in these systems 
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will allow us to correlate the observed properties with the degree of orbital alignment, and thus 
to search for the physical mechanisms causing observed misalignments. We will also have the 
information to be able to begin to back-out the physical angle between the spin and orbital angular 
momentum vectors, not just the projection of these angles on the plane of the sky. 

Funding for this Discovery mission is provided by NASA's Science Mission Directorate. Many 
people have contributed to the success of the Kepler Mission, and it is impossible to acknowledge 
them all. Valuable advice and assistance were provided by Willie Torres, Riley Duren. M. Crane, 
and D. Ciardi. Special Technical help was provided by Carly Chubak, G. Mandushev, and Josh 
Winn. We thank E. Bachtel and his team at Ball Aerospace for their work on the Kepler pho- 
tometer and R. Thompson for key contributions to engineering; and C. Botosh, and J. Fanson, 
for able management. GWM thanks acknowledges support from NASA Cooperative Agreement 
NNX06AH52G. 
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Table 1. System Parameters for Kepler-8 



Parameter Value Notes 



Transit and orbital parameters 

Orbital period P (d) 3.52254^°;°°°^ A 

Midtransit time E (HJD) 2454954.11821° °°°^ A 

Scaled semimajor axis a/ij* 6.97^q'24 A 

Scaled planet radius R P /R+ - 09809 io.ooo46 A 

Impact parameter b = a cos i/R* 0.724 ± 0.020 A 

Orbital inclination i (deg) 84.07 ± +0.33 A 

Orbital semi-amplitude K (ms^ 1 ) 68.4 ± 12.0 A,B 

Orbital eccentricity e (adopted) A,B 

Observed stellar parameters 

Effective temperature T cff (K) 6213 ± 150 C 

Spectroscopic gravity log g (cgs) 4.28 ±0.10 C 

Mctallicity [Fe/H] -0.055 ± 0.03 C 

Projected rotation Velocity v sin i (kms- 1 ) 10.5 ±0.7 C 

Absolute (Helio) radial velocity (kms" 1 ) -52.72 ± 0.10 B 

Derived stellar parameters 

Mass M*(M ) 1.213±g;g|| C,D 

Radius R±(Rq) 1.486j:°;°g C,D 

Surface gravity log (cgs) 4.174 ±0.026 C,D 

Luminosity L* (L Q ) 4.03±°,;54 C,D 

Absolute V magnitude M v (mag) 3.28 ± 0.15 D 

Age (Gyr) 3.84 ±1.5 C,D 

Distance (pc) 1330 ± 180 D 

Planetary parameters 

Mass M P (Afj) 0.6031°;^ A,B,C,D 

Radius R P (Rj) 1.4191° °^ A,BC,D 

Density p P (gem" 3 ) 0.261 ± 0.071 A,B,C,D 

Surface gravity log g P (cgs) 2.871 ±0.119 A,B,C,D 

Orbital semimajor axis a (AU) °- 0483 io:ooi2 E 

Equilibrium temperature T cq (K) 1764 ± 200 F 



Note. — 
A: Based on the photometry. 
B: Based on the radial velocities. 

C: Based on an SME analysis of the Keck-HIRES spectra, 

D: Based on the Yale-Yonsci stellar evolution tracks. 

E: Based on Newton's version of Kepler's Third Law and total mass. 

F: Assumes Bond albedo = 0.1 and complete redistribution. 
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Table 2. Relative Doppler Velocity Measurements of Kepler-8 



HJD 


RV 


(TRY 


(-2450000) 


(ms- 1 ) 


(ms- 1 ) 


4984.040 


-22.15 


34.5 


4985.043 


52.45 


34.5 


4986.065 


-42.39 


34.0 


4987.060 


-74.37 


33.6 


4988.027 


63.04 


33.7 


4988.973 


98.87 


33.8 


4995.103 


46.58 


34.2 


5014.891 


-64.61 


33.7 


5015.975 


-29.97 


33.4 


5017.020 


-28.28 


33.3 


5075.788 


41.65 


33.9 


5109.850 


-122.00 


32.8 


5110.790 


-67.42 


31.9 


5133.839 


49.58 


32.8 


5133.718 


58.11 


32.6 


5133.732 


51.13 


32.5 


5133.747 


71.85 


32.7 


5133.761 


34.51 


32.8 


5133.776 


24.38 


32.6 


5133.790 


-9.49 


33.2 


5133.805 


-31.04 


33.2 


5133.821 


-14.32 


32.6 


5134.740 


-23.48 


31.9 


5135.790 


48.50 


32.6 
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Fig. 1. — The detrended light curve for Kepler-8. The time series for the entire Kepler photometric 
data set is plotted in the upper panel. The lower panel shows the photometry folded by the period 
P = 3.52254 days, and the model fit to the primary transit is plotted as a solid line. Data for even 
and for odd transits are denoted by '+'s and '*'s, respectively. The even and odd transits do not 
exhibit significant variations in transit depth indicative of the primary and secondary eclipses of an 
eclipsing binary. The photometry at the predicted time of occultation of the planet is shown just 
above the phased transit curve, showing no significant dimming at the level of 0.1 mmag. The lack 
of any dimming at the predicted time of occulation, along with the absence of both astrometric 
motion in and out of transit and the absence of detected stellar companions, leaves the planetary 
interpretation as the only plausible model, confirmed by the Rossiter-McLaughlin detection. 
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Fig. 2. — The image of the field near Kepler-8 obtained from the Keck 1 telescope guider camera 
with a bg38 filter, making the bandpass approximately V and R combined, similar to that of the 
Kepler photometer. North is up, and east to the left, and the plate scale is 0.397 "per pixel. The 
main star, Kepler-8 is the brightest source at the center of the image. A faint star is located 
3.8 "to the northwest of Kepler-8, having a flux 0.75% that of Kepler-8 and residing within the 
Kepler photometric aperture. With a flux less than 1% that of Kepler-8, the companion cannot be 
responsible for the periodici 1% dimming, even if it were a background eclipsing binary. But this 
neighboring star is sufficiently bright to account for the astrometric displacement of the Kepler-8 
that occurs during transit of its planet, explaining the shift of the photocenter of the two stars by 
0.0001 Kepler pixels. 
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Fig. 3. — Flux(upper) and astrometric centroid (middle, lower) timeseries of Kepler-8b from Kepler 
quarter-1 data. A high-pass filter has been applied to each time series to remove non-transit 
signatures on timescales longer than 2 days. There is no obvious displacement of the centroid of 
light in either the row or column residuals. This lack of motion indicates that any background 
eclipsing binary would need to be located within ~0.001/.01 = 0.1 pixels corresponding to 0.4 "of 
the target star in order to explain the photometric transit signals but avoid showing motion. 
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Fig. 4. — Flux versus residual astrometric centroids in row (x) and column (+) for Kepler-8. The 
measurements are the same as those presented in Figure El This "rainplot" shows no evidence that 
the centroids are systematically shifting during each transit event, which would be indicative of a 
blend with a background star. If this were the case, identification of the actual star responsible for 
the transit features in the light curve would require inspection of a high resolution image or catalog 
to unravel the mystery. Since the in-transit points "rain" straight down, rather than slanting along 
a diagonal, any blend scenario would require a very small separation on the sky between the two 
sources. 
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Fig. 5. — The results of a "poor man's" difference image analysis are presented in this figure: (A) 
(upper left) is the result of subtracting the average of the in-transit images from the average of a 
similar number of just out-of-transit images for one of Kepler-8b's transits. The pixels summed to 
form the photometric brightness measurements are indicated by the blue circles. For comparison, 
the average of the out-of-transit images are displayed in (B - upper right) and (C - lower left). The 
raw flux light curve from QO is displayed panel (D - lower right). The difference image, Panel (A) 
is remarkably similar to the direct images. Had the transits been due to a background eclipsing 
binary that were 0.5 or more pixels offset from the target star, then the difference image would 
display a noticeable shift in the location of the pixel containing the peak energy. Instead, both 
difference and direct images supporting a coincident localization of the source of the transits to the 
target star, Kepler-8. 
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Fig. 6. — The orbital solution for Kepler-8. The observed radial velocities obtained with HIRES on 
the Keck 1 Telescope are plotted together with the best-fit velocity curve for a circular orbit with 
the period and time of transit fixed by the photometric ephemeris. Uncertainties are large due to 
the faintness of the star (V=13.9 mag), the high rotational Doppler broadening (v suu= 10.5 ±0.7), 
and contamination from background moonlight in most spectra. The measured velocities vary high 
and low as predicted by the orbit from the photometric transit curve, supporting the planet model. 
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Fig. 7. — The Rossiter-McLaughlin effect for Kepler-8. The observed radial velocities obtained 
with HIRES on the Keck 1 Telescope are plotted together with the predicted velocities from the 
best-fit model of the Rossiter-McLaughlin effect. The high velocities near ingress and low velocities 
near egress imply a prograde orbit. The asymmetry in the model curve indicates a projected tilt 
of the orbit relative to the star's equator, as well as a transit with a large impact parameter on the 
star. 
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Fig. 8. — A sketch of the transit geometry for Kepler-8. The Kepler photometry and R-M 
measurements set the impact parameter, b = 0.724 ± 0.020, and the angle projected on the sky 
between the star's spin axis and the normal to the orbital plane, A = —27°, indicating a prograde 
orbit with a moderate inclination. 



